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CH ' Abstract. 

In this contribution, we discuss three situations in which complete 
integrability of a three dimensional classical system and its quantum version can 
be achieved under some conditions. The former is a system with axial symmetry. 
In the second, we discuss a three dimensional system without spatial symmetry 
which admits separation of variables if we use ellipsoidal coordinates. In both 
cases, and as a condition for integrability, certain conditions arise in the integrals 
C^ i of motion. Finally, we study integrability in the three dimensional sphere and a 

particular case associated with the Kepler problem in S 3 . 



1. Introduction. 

in 



For a system whose configuration space has dimension n, integrability in the sense of 
Liouville- Arnold requires the existence of a number n of constants of motion which 
should be in involution (the Hamiltonian itself is one of them). Integrable systems 
are rare, yet many systems which are important from a physical standpoint turn out 
r — ' to be integrable, with the outstanding examples of the harmonic oscillator and the 

Kepler problem, which in fact are superintegrable (they have more than n functionally 
independent constants of motion, albeit not all of them are in involution). 

A type of integrability appears to be particularly relevant: the quadratic 
integrability. A standard natural Hamiltonian has a kinetic part which is quadratic 
in the momenta. If all n — 1 additional constants of motion have a similar kind of 
structure, we will speak of quadratic integrability (understanding quadratic as 'at 
most quadratic', i.e., allowing possibly for constants of motion which are linear in the 
momenta). 

For a 2D system in any constant curvature configuration space, quadratic 
integrability means the existence of a single constant of motion quadratic in the 
momenta; this case is simpler as there is no any extra condition ensuring the 
additional constants are in involution [Tl 126] . For any space of constant curvature, 
integrable systems have a Hamiltonian which is separable in confocal (or general 
elliptic) coordinates in that space [2] [5] . These coordinates include generic ones (as 
the elliptic coordinates in the sphere S 2 or the Euclidean plane E 2 ) as well as all its 
possible degenerations or limiting cases. This also holds in the 3D case [21 |S] ■ It is 
generally believed that a similar connection among quadratic integrability and the nD 
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version of confocal coordinates does somehow hold, but to the best of our knowledge 
a proof is not available in the literature. 

In this work we explore some aspects of quadratic integrability for a system in 3D 
configuration spaces E 3 or S 3 . The results should be considered as a stage towards 
studying how precisely integrable systems in Euclidean or in any constant curvature 3D 
space are related to separable systems in the generic version of ellipsoidal coordinates, 
associated to confocal quadrics [TBI HZ]. We discuss mainly the Euclidean case, but 
extension of these results to spaces with any constant curvature seems to be possible, 
and as a hint in that direction, we discuss the integrability with axial symmetry and 
the Kepler problem in a space with positive constant curvature S 3 (see also [13]). 

2. 3D systems with axial symmetry. 

In this first part, we consider a three dimensional system in Euclidean space with 
axial symmetry around an axis given by the unit vector n = (m, n 2 , nz). The classical 
analysis of this system makes use of canonical coordinates (p, x) = (pi,Xi), i = 1, 2, 3 
and Poisson brackets iA,B\ = -S^-4^- ~ 4^-4^-- We want to characterize those 

l ' ' ap k dx k dx k dp k 

systems with axial symmetry which are completely integrable. For these systems 
dynamics is described by a Hamiltonian of the form 

tf=f^(x), (1) 

where U(pc) is a time independent potential. Axial symmetry requires that the angular 
momentum L around the symmetry axis 

L = n • (x x p) (2) 

be a constant of motion, i.e., {H, L} = 0. Complete quadratic integrability requires 
the existence of a third independent constant of motion Hi, which is in involution 
with H and L, i.e., {H, Hi} — {L,H\\ — and has a 'Hamiltonian form' in general 
position dependent yet quadratic in the momenta, plus a time independent 'potential', 
with form 

Hi = ^ Pi g ik (x)p k + $(x). (3) 

The tensor g lfc (x), which determines the 'kinetic' part of the constant of motion Hi 
has to be a Killing tensor for the Euclidean metric. Its determination can be done 
by assuming that the 'kinetic' terms in {]]) and ([3|) commute among themselves and 
also with the angular momentum L, with respect to the Poisson bracket. If we choose 
the angular momentum direction as the z axis, then the most general expression for 
Vl g lk (x) Pk is [20]: 

Pl g tk (x) p k =L\ + Ll + apj + /3(piL 2 - p 2 Li) + jp 3 L 3 (4) 

to which of course any linear combination of the quadratic Casimirs of the Euclidean 
algebra, p\ + p\ + p\ and piLi + p 2 L 2 + P3L 3 can be added. Here we will discuss 
mainly the family with a < 0, /3 — 0, 7 — 0, for which 

Pi g ik (x)p k = L 2 i + Ll-a 2 pl (5) 

with a a positive constant. If now the angular momentum direction is along a general 
unit vector n = (711, n 2l n 3 ), by direct computation we get for the Killing tensor g lfc (x) 
the expression 

g ik (x.) = <f fc (x • n) 2 - x • n( Xl n k + x k n t ) + (x 2 - a 2 )n t n k , (6) 



Classical and quantum integrability in 3D systems. 3 

For more details, see |5]. The commutation of L with H and H\ restricts the form of 
the functions U and $, collectively called 'potential' terms to: 

[/(x) = [/(x 2 ,(x.n) 2 ), $(x) = $(x 2 ,(x.n) 2 ) (7) 

(note that the condition of axial symmetry around the n axis implies that U and $ 
cannot depend on the azimuthal angle <p) and the commutation of H and H\ leads to 
the equations: 

d i 3>(x.)=g ik (x)d k U(x), i J * = l,2,3, (8) 

where we always assume summation over repeated indices. In order to obtain the 
solutions of ©, we diagonalize the matrix with components g tk (x). Its eigenvalues 
A(x) and eigenvectors A(x) are obtained from the matrix equation 

(g ik (x)-\(x)S ik )A k (x) = 0. (9) 

The eigenvalue equation © has the following solutions: 

x 2 




A±(x) = — ^— ± W I — ^— j + a 2 xn 2 ; A (x)=xn^,(10) 

with corresponding eigenvectors given by 

Afw = a,A T (x) , a°(x) = x "_ x x x n2 = a^(x) , (ii) 

where the index i = 1,2,3 labels the component of the corresponding vector field and 
ip is the azimuthal angle around the n-axis. These vectors form a new coordinate basis, 
{VA_, VA + , Vip} referred to which the matrix g zk (x) is diagonal. Taking into account 
that, due to the geometric symmetry, U and $ cannot depend on tp, JHJ) becomes 

d,x+ 9+$ + a,A_ a_$ = g lk (d k x+ d+u + a fe A_ d_u) , (12) 

where d± stand for or - • In the basis {VA_, VA + , V<^} that diagonalizes the matrix 
g lfc (x), equation (fT2"|) decouples in 

9+$ = A_a+C/, 9_$ = X-d+U. (13) 

This means that 

$-A+t/ = -/(A+), $-A_C/ = -. 9 (A_), (14) 

where /(A+) and <?(A_) are arbitrary functions on their respective variables. These 
equations give an expression for the 'potentials' as follows: 

u.iM^L, *^ /(A+) -^- 9( ,_,, 15 , 

These are the most general expressions for the potentials f (x) and $(x) compatible 
with our choice for H i and with the requirement of H , _ff i , L being in involution. 

A comment is here in order: the presence of a constant of motion which is first- 
order in the momenta, i.e. is a Noether constant, means that the system is invariant 
under rotations around the axis n. Then this situation reduces to one that is essentially 
two-dimensional, with a single azimutal coordinate tp added to the ordinary 2D elliptic 
coordinates in a fixed plane containing the axis n. In the 3D Euclidean space these 
coordinates are the spheroidal (oblate) cordinates [55], which play the role of 3D 
separation coordinates for the system, as we will show explicitly in the next section. 
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2.1. Integration. 

The integration of the classical system with axial symmetry with respect to the n 
axis is based in the following idea: If e, = h(p,x), i = 1,2,3 are three independent 
functions with {ei,ej} — 0, then, there exists a function -F(x, e) with e = (ei,e2,e3) 
and 

3F(x,e) 

Pk = n ■ 16 

ox k 

The function -F(x, e) is the characteristic function in the Hamilton-Jacobi approach. 
In our case, the chain rule and (fT6| gives the following expression for the 
momentum p ~ (pi,P2,P3)'- 

p = VA+ d+F(X+, X_,tp) + VA_ 9_F(A+, A_, tp) + 2 " XX , 2 d v F(X + , X-,tp) .(17) 

x^ — (x ■ n) z 

Then, the kinetic part of H can be written in the following form: 

p 2 = (VA+) 2 (d+F) 2 + (VA_) 2 (cLF) 2 + -5— £ = 2m (E-U), (18) 

x- 5 — (x • n) z 

where I is the value of the integral of motion L corresponding to the angular 
momentum around n, i.e., 

(n x x) ■p = n-L = l = d ip F(X+,X-,tp) (19) 

and E is a given value of the constant of motion H . 

For the kinetic term of Hi, we obtain a similar expression: 

Pl g lk Pk = A_(VA + ) 2 (cV^) 2 + A + (VA_) 2 (5_F) 2 + 2 X, " )2 ^, 2 = 2m(E 1 - *) , (20) 

x^ — (x • ny 

where E\ is a constant value of FL\ . 

With some calculations [9], we obtain the partial derivatives of the function 
F(X + , A_, tp) in terms of its arguments: 

-.2 1 L „ fr , rrt ,„ „„ A+Z 2 



(d+FY = - 2m(\ + (E -U)- (E 1 - $)) - —^ , (21) 

4A+(A + + a 2 ) ^ A + + a- ' 



< 3 - f »' - 4>_(>- +«') { 2, " (A - (E " "' " (E ' " *» " A^} - <22> 

Due to (|14[) . pT|) does not depend on A_ and (I22p does not depend on A+. From 
this simple idea, we conclude that the function .F(A+, A_, tp) is of the form: 

F(X+,X-,tp)=A(X + )+B(X_) + etp, (23) 

where A(-) and B{-) are functions of one variable only. Thus, we have separated 
variables in the Hamilton-Jacobi generating function. This permits us to obtain the 
equations of motion in terms of the variables (A+, X—,tp). The final result is 9J 

A, A_ 4 



X+d + A(X+) A_d_5(A_) m 



+ , „ ,, n ,„,, ; =0 



A + (A + + a 2 )<9 +J 4(A+) A_d_(A_+a 2 )B(A_) 



Za 2 / A+ A_ 



4 U + (A + +« 2 ) 2 M(A+) A-<9-(A-+a 2 ) 2 B(A-) 



= 0(24) 
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where the upper dot represents the derivative with respect time. These equations are 
at least formally integrable. 

2.2. Quantum, case. 

Canonical quantization of the functions H,L,Hi give respective operators that we 
represent with the same symbols. Then, complete integrability means that 

[H,L] = [H 1 ,L] = [H,H 1 ]=Q, (25) 

where [A, B] — AB — BA denotes the commutator of the operators A and B. Equation 
(|25| implies the existence of simultaneous eigenfunctions ip(x) of these three operators: 



(H - E)4' = (^ - E^j = {L- £}4j = , (26) 
or equivalently, 

- AV>(x) = 2m{E - U(x))ip(x) (27) 

-A 1 V'(x) = 2m{E 1 - $(x))V»(x) (28) 

-id (p ip(x)=ti>(x), (29) 

where, 

A = d k d k , A 1 =d j g jk (x)d k . (30) 

We can express (f27|) and (j28|) in terms of the variables A±. Using (|29|) we obtain 
respectively, 



4A + (A++a 2 )V++ + 2(a 2 +3A + )^+ - -*-^ </> 



l 2 X. 



2m[(E-U)\+-(E 1 -$)]Tp{3l) 



4A_(A_+a 2 )?A— + 2(a 2 +3A_)V>- - t ^ 

A_+cr 



2m[(E-U)X--(Ei-<S>)]ip .(32) 

Note that (fBTj) and (|32[) depend only on A + and A_ respectively. Then, the wave 
function ip(x) can be factorized as 

V>(x)=V + (A + )VT(A-) e ^ (33) 

so that (|3"0)) , (|3"Tj) and (f2"5)) can be written as equations depending solely on the variables 
(A+, A_, ip) respectively. Once we have solved equations (|30|) and (|3Tj) . we have solved 
the problem of finding solutions of (|26|) or, equivalently, (|27H29|) . In order to finding 
solutions, we have to give explicit expressions for the potentials U and $. These 
expressions must be obtained from (|15|) and the form of the functions /(A+) and 
S(A-). 

In order to illustrate the search for solutions of (|3"Tj) and ([32]) . let us propose a 
simple although nontrivial choice of /(A+) and g(X-) so that ([30)) and (f3Tj) are exactly 
solvable. Here, we propose 

/(A+) = $ - X+U = , 5 (A_) = $ - A_C7 = -Q(A_ + a 2 ) , (34) 

where Q is a constant. This choice give the following form for the potentials: 

tf(x) = - g£±£, *(x) = -Q^-±^. (35) 
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The steps to solve (15U|) and (|3ip with these potentials are the following: i) From 
(fTU)l. we conclude that A + is always positive meanwhile the values of A_ lie on the 
interval [— a 2 ,0]. This suggest the following change of variables: A + = a 2 sinh 2 a and 
A_ = —a 2 sin f3. The new coordinates (a,f3,<p) are known as the oblate spherical 
coordinates. This change of variables transforms (|30|) and (|3Tj) into two new wave 
functions |9]. ii) The new change of variables t := sinh a and u := sin/3 and the 
introduction of the new parameters, £ := 2ma 2 E, q := 2ma 2 Q, £\ :— £ 2 + 2mEi, G := 
£\ + £ and q' = q — £, with t = ia transform equations (|30p and (|3ip respectively into 

,„ 9x d 2 ib + (a) „ d%b + (a) ( I 2 „, 9 , ] ,,. , ,„„. 

(!-» )^3- J - 2 «^Z LJ + G -T— 3- f (!-" ) K» ' ( 36 ) 



da 2 da 

2 d?il>-(u) di/j-(u) f I 2 ,, 2 ,\ , , ,„_, 

(1 - u) ^i 2 -- 2u ^^ + i G -T-^ + ' ?(1 - u) r (u)= (37) 

Equations (|36ti37|) are spheroidal equations, a type of second order differential 
equation which has been studied [5T] . A brief discussion of the solutions for (|36H37[) 
and therefore for (|3HI32p can be found in [}J]. 

3. General 3D integrable systems. 

In this section, we intend to generalize the previous discussion to the case in which no 
(Noether) symmetry is present. As in the previous case, we shall discuss the classical 
point of view first and then its quantum counterpart. The key of the solution to this 
general case will be the choice of a proper system of coordinates, as we shall see. This 
justifies the next subsection. 

3.1. The ellipsoidal coordinates. 

The generic ellipsoidal coordinate system in the Euclidean three dimensional space is 
determined as follows. Fix three positive numbers a, b and c with a > b > c, and 
consider the one-parameter family of quadrics with equation 

2 2 2 

-1, (38) 



a 2 + £ b 2 +£ c 2 +i 

Through any given point (x,y,z) in Euclidean 3D space, there passes precisely three 
such quadrics, corresponding to the values A,/i, v of the parameter £ which lie in the 
intervals 

- a 2 < v < -b 2 < (j, < -c 2 < A . (39) 

One quadric is an ellipsoid (because a 2 + A > b 2 + A > c 2 + A > 0) 

2 2 2 

X V Z 

a + A cr + A <r + A 
other is an one sheeted hyperboloid (because a 2 + /i > b 2 + \x > > c 2 + /j) 

2 2 2 

x y Z 

—o + 15* + 3 = 1) ( 41 ) 

a z + fj, b z + /! c z + /! 

and the last is a two sheeted hyperboloid (because a 2 + v > > b 2 + v > c? + v) 

X " • ^ • ^ =1. (42) 



b 2 + v c? + v 
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Solving (|40j|42p for x 2 , y 2 and z 2 , we obtain the parametrization of 3D Euclidean 
space in terms of ellipsoidal coordinates, A, /i and v. 

2 _ (A + a 2 )( M + a 2 ){v + a 2 ) , , 

(b 2 - a 2 )(c 2 - a 2 ) ' [ } 

, _ {\ + b 2 )^ + b 2 )(v + b 2 ) 

V (~1 V>.\( „2 J.2N ' y^> 



(c 2 - b 2 ){a 2 - b 2 ) 
2 _{\ + c 2 ){i, + c 2 ){u + c 2 ) 



(45) 



(a 2 - c 2 )(b 2 - c 2 ) ' 

Inversion of (|43M45[) gives ellipsoidal coordinates in terms of Cartesian coordinates. 
Note that ellipsoidal coordinates are only one-to-one in the interior of each octant of 
the Euclidean space determined by the three coordinate planes x = 0, y = 0, z = 
through the origin. 

3.2. Classical 3D integrable systems. 

Let us study the situation from classical point of view first. Consider a three 

2 

dimensional system with Hamiltonian given by H = #— + £/(x) and look for two 
additional independent integrals of motion of the form: 

Hi = ^ Pi g\ j (x) Pj - + C/i (x) , H 2 = ^- Pi gl j (x) Pj + U 2 (x) , (46) 

with {H, Hi} = {H,H2} — {H\,H^] = 0, where again the brackets are the Poisson 
brackets. The tensors g\ J (x) and g% (x) can be obtained using the hypothesis of 
commutativity of the kinetic parts of H, H\ and H2 (due to their structure, we could 
extend the name 'Hamiltonian' to Hi, i = 1,2 too). The condition 

{ Pl gl k (x)p k ,p 2 } = (47) 

implies that Pig\ k (x)p^ should be a quadratic function of the momenta. A particular 
choice follows from a further restriction of Hi , Hi being invariant under any reflection, 
so that both can be expressed as a linear combination of the squares L 2 of the 
components Li of the angular momentum and p 2 of the linear the momentum pi. 
There is some arbitrariness in the choice of Hi , for which we set 

p l9 \ k {*)p k := L 2 -p\(b 2 + c 2 )-p 2 (a 2 + c 2 ) - p 2 (a 2 + b 2 ) , (48) 

where a, b and c are the constants which specify the system of ellipsoidal coordinates, 
see (|39l) . In matrix form, the tensor g\ k (x) defined by (|4"5)l can be written as 

(z 2 + y 2 — (a 2 + c 2 ) -xy —xz \ 

-xy x 2 + z 2 - (a 2 + c 2 ) -yz . (49) 

— xz —yz x 2 + y 2 — (a 2 + b 2 ) J 

Analogously, the commutativity of the kinetic part of Hi with those of H and 
Hi, 

{P % gf{*)p k , Pl gl k (x)p k } = (50) 

fixes the linear combination of the components of L 2 and of p 2 in Hi as 

Pi gi k (x) Pk = -(L 2 a 2 + L 2 b 2 + L 2 c 2 ) +p 2 b 2 c 2 +p 2 a 2 c 2 +p 2 a 2 b 2 , (51) 



with matrix form given by 




/ -c?y 2 - b 2 z 2 + b 2 c 2 


c 2 xy 


gf(x) = c 2 xy 


—<?x 2 — a 2 z 2 + a 2 c 2 


\ b 2 xz 


a 2 yz 
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b 2 xz 

a 2 yz \ .(.-,2) 

-b 2 x 2 - a 2 y 2 + a 2 b 2 

Matrices P§|) and (|5^|) are symmetric and therefore diagonalizable. In addition, 
these matrices commute and therefore admit a diagonal form in the same basis. This 
basis is given by 

E x := (d x X, d y \, d z X) , E^ := (d x fi, d y /j,, d z y) , E u := (d x v, d y v, d z v) . (53) 

These vectors are mutually orthogonal [10]. Note that A, fj, and v are functions of 



the Cartesian coordinates, so that the partial derivatives in (|53[) make sense. The 
eigenvalues of (|49[) and (|5^|) are easy to obtain and are given in the following table: 

ff «(x) gl j (x) gi j (x) 

Ex 1 n + v \xv (54) 

Eft 1 X + u \v 

E v 1 A + fi Xfi 

Thus, we have determined the kinetic parts of H\,Hi in terms of ellipsoidal 
coordinates. Our next goal is to obtain the most general form of the 'potential' terms 
f(x), £/i(x) and [^(x). From the commutation relations {H, Hi} = {H, H2] = 
{Hi, H2} = 0, we obtain the relations, 

d*Ui=g?d k U, d*U 2 =gfd k U, gfd k U 2 = gfd k Ui . (55) 

Then, after a calculation that makes use of the chain rule involving ellipsoidal 
and Cartesian coordinates, the form (|53|) of the eigenvectors of matrices g\ k (yz) and 
g 2 k (x) and the orthogonality of these vectors, we can obtain the most general form of 
the potentials U, Ui and U 2 satisfying {H,Hi} = {H,H 2 } = {H X ,H 2 } = [10 . This 
is: 

U = IW + MP) + <y) ^ 

(\-u)(X-fj,) (fi - v)(jj, - X) (v - fi)(u - X) ' 

(,u + v)l(X) , (v + X)m(n) (X + (jL)n(v) 

1 (A-^)(A- Ai ) + (^-z,)0i-A) + (v-Mv-X)' ( > 

u = (V)IW + {vX)m{n) + (Xn)n(v) 

(X-u)(X-n) (ji-v)(n-X) {v - n)(v - X) ' 
where l(X), m{ji) and n(v) are arbitrary functions of their arguments. This is the 
most general form of the potentials compatible with complete integrability. 

3.3. The quantum case. 

We shall briefly comment the procedure here; details can be found in [TU]. First of 
all, we obtain the quantum operators H, Hi and H 2 by canonical quantization of 
their classical counterparts (l46l) . Then, complete integrability means that these three 
Hamiltonians commute with each other. Therefore, we can find wave functions "0( x ) 
such that 

{H - E)if> = (Hi - Ei)i/t = {H 2 - E 2 )i/j = . (59) 



2 + 2m[X 2 {t)E - \{t)E 1 +E 2 ~ Z(A(i))]V = , (62) 
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We can prove that the use of ellipsoidal coordinates A, [i and v, the form for 
the potentials given in ([55|) and the factorization t/>(x) = ijj(\)(t)(n)Lp(v) show that 
equations (|59|) are equivalent to three wave functions solely in the variables A, [i and v 
respectively, so that complete separation of variables is also achieved in the quantum 
case [10]. The equation for A gives (ijj\ denotes the derivative of ij) with respect to A): 

4(A + a 2 )(A + 6 2 )(A + c 2 )V^aa + 2[a 2 6 2 + a 2 c 2 + b 2 c 2 + 2\{a 2 + b 2 + c 2 ) + 3A 2 ]^a 

+2m(X 2 E- XE 1 + E 2 - l(X))tp = . (60) 

The other two equations are just obtained by replacing A by ^ and v respectively. 
Equation (|60p can be written in a more compact form by means of an of a change of 
variable to the auxiliary variable t. This is given by means of the condition: 

A'(i) = 2v/(A + a 2 )(A + & 2 )(A + c 2 ). (61) 

In terms of t, equation ([HO)) has the follwing form 

~dl 

where 1(A) has already appeared in (|56ti58|) . In addition to (|62|) . there arc two 
other equations, one for /x and the other for v. These three equations are similar 
and, in particular, all depend on the function X(t) even in the free particle case, 
/(A) = m(n) = n{u) = 0. 

4. Integrability in spaces of constant curvature. 

There are three types of homogeneous three dimensional manifolds with constant 
curvature and Riemannian positive definite metric, which are the sphere S , the 
Euclidean plane E 3 and the hyperboloid H 3 (for general references on integrability in 
spaces of constant curvature, see [5Dj). In this section, we shall deal with the standard 
S 3 with curvature equal to 1 . As is well known, S 3 can be realized as the submanifold 
X 2 + X 2 + X 2 + X 2 = 1 in a 4D Euclidean ambient space, with the induced metric. 
A convenient coordinate system for S 3 is provided through stereographic coordinates, 
the point on S 3 parametrized by x e R 3 U oo being given by: 

X = (X,X 4 ) = (-^, %~4)- (63) 

\x J +x * +X J 

where x is a parameter which should be different from zero. With this choice x = 
corresponds to the sphere's 'South Pole' (0, 0, 0, —1), and the projection is made from 
the North Pole (0, 0, 0, 1) to the plane X4 = 1 — \- No actual generality would be lost if 
we fix a particular value for %; the preferred choice \ = 2 corresponds to projecting over 
the plane tangent to the sphere's south pole and expressions are generally more clear 
for this choice, because in this case near the South Pole the coordinates x approach the 
ordinary Cartesian coordinates in the tangent E 3 with the same scaling as lengths on 
the sphere (a rescaling is required for other values of x) and neglecting terms which are 
higher order in x brings the Euclidean expression with the correct factors, as obvious 
for instance in (l65 |) -([66 ]) . 

The Lagrangian for the free evolution on S 3 can be taken in terms of the ambient 
space coordinates as: 



L Q = 



y (* 2 + if) , (64) 



Classical and quantum integrability in 3D systems. 10 

which can be expressed in stereographic coordinates as 



Pi = m 7^ t rsw ^- ( 66 ) 



m x 2 . 

£ = 7 4 X 2 - 65 

Alternatively, we may start from the angular momentum tensor in ambient space, 
with components M a p :— X a Xp — X a Xp and define the free lagrangian in S 3 as 
(proportional to) the square of this angular momentum tensor; this leads to the same 
free lagrangian. The canonical conjugate momenta associated to the stereographic 
coordinates Xi is 

' (x 2 + X 2 
and the Legendre transformation provides the following free Hamiltonian: 

The symmetry group for S 3 is SO (A) and (|67[) is invariant with respect to the 
sphere isometries in 50(4). Then, H has 6 'kinematic' integrals of motion, which are 
generators of 5*0(4); in terms of the stereographic parametrization these generators 
are: 

Li = eijkXjPk , Ki= — \2xiTp ■ x - pi(x 2 - x 2 )) , (68) 

with Poisson brackets 

{Li, Lj} = -eij k L k , {K t , Kj} = —et-jkLk , {Li, Kj} = -6ij k K k .(69) 

(these close also an so(4) algebra as expected). There are two Casimirs given by 
K 2 + L 2 and K • L. The free Hamiltonian is proportional to the first Casimir: 

ff = ^-(K 2 +L 2 ). (70) 

while the second Casimir vanishes. 

Then, we are going to pose in S 3 the same question studied for the flat case in 
section 1. Assume the existence of a symmetry axis, that, for the moment, we can 
identify with the z-axis. Then, consider a Hamiltonian with kinetic term given by Hq 
plus a time independent potential U(x), 

ff=-^(K 2 +L 2 ) + C/(x). (71) 

Zm 

Then, if this system has to be completely integrable, we need finding a new 

'Hamiltonian' H\ with 'kinetic' part T\ = (2m)~ 1 pig l: > (x)pj and 'potential' ^(x), 

Hi=T + $(x), independent of H and such that {H, L 3 } = {H x , L 3 } = {H, H x ] = 0. 

As in Section 1, the kinetic term T\ is obtained under the condition that it commutes 

with Hq and L 3 . There is an S 3 analog of the most general 'kinetic' term $$, but we 

simply write down the version for S 3 of the Euclidean Killing tensor ([5]) leads to the 

constant of motion: 

Hi = ^ (Lj + Lj- a 2 Kl) + $(x) , (72) 
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where a is a constant playing a role fully similar to those in (J5]) . If the symmetry axis 
were arbitrary in the direction given by the unitary axis n = (ni, n 2 , n 3 ), the matrix 
g lJ has the following form in terms of the stereographic parametrization: 

g lJ {~x) = Sij(x ■ n) — x • n(xinj + Xjiii) + UiUjX 

2 

- ^(2x,x • n - n ( (x 2 - X 2 )) (2x 3 x • n - n,(x 2 - X 2 )) -(73) 

The requirement {H, Hi} = leads to the version for S 3 of equation IJ5J), which 
in this case reads 



9,$(x)=G^(x)5,[/(x), with GV(x)= ,J ,\; . (74) 



g ij "(x) 

(x 2 + x 2 



We shall denote the eigenvalues of G y (x) by {/i+, /i-,/i3}. The two former are 
the roots of the quadratic equation z 2 — Pz — Q = 0, where 

_ x 2 - Q 2 [(x 2 - X 2 ) 2 +4(x.n) 2 ] _ (x.n) 2 a 2 

(x 2 + X 2 ) 2 ' (x 2 + l) 2 l j 

and /i3 = Q/a 2 . As in the flat case described in section 1, the eigenvalues /j,± together 
with the azimuthal angle (p define a coordinate system, which is the analogous for 
the sphere S 3 of the spheroidal oblate coordinates [23]. As in the flat case, their 
corresponding eigenvectors have components di/i- and 9i/i+, i,j = 1, 2, 3 respectively, 
i.e., 

G lJ dji.i- = [i+difi- , G lJ dj^ + = pt,-diH + . (76) 

Finally, following a procedure similar to that studied in section 1, we can show 
that equations (fTB")) provide the most general form for the 'potentials' U(x) and $(x) 
on the variables fi±. The final result is: 

jj = /(M+) -g(M-) $ = M-/(^+)-M+g(^-) / 7? n 

M+ - M- M+ - M- 

where f((x+) and g{fi-) are arbitrary functions. Note that symmetry prevents that 
f and $ to depend on the azimuthal angle <p, as constancy of L% means invariancc 
under rotations with z-axis. 

4--1- Application to the Kepler problem. 

In this subsection, we shall discuss a particular case of very special importance: 
the Kepler problem in the standard sphere S 3 with curvature equal to 1, which 
provides an example of the type of situation just described (recall a factor which 
dimensionally is a lenght square has been taken as equal to 1, and hence is invisible 
in the expressions). For more details on the Kepler problem in spaces with constant 
curvature, see H| [ED US E3 [25] . 

In the sphere S 3 (and also in any dimension n > 3) , this system is an example of 
a maximally superintegrable system. Restricting attention to the 3D case, the Kepler 
Hamiltonian is 

H = J- ( K 2 +L 2 )--4t (78) 

2m tan(r) 

where the kinetic part is (|7ip and the potential term depends only on the intrinsic 
distance r in S 3 to the potential center and k > for the atractive case. The potential, 
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with center at sphere's south pole, can be expressed in terms of the ambient space 
coordinates as 

''' =k , Xi (79) 



tan(r) j X \ + X* + Xf 

and in terms of stereographic coordinates the Kepler Hamiltonian is: 

ff= J_ (K * + LV^^. (80) 

where the generators in the kinetic part are given in (|68[) . This Hamiltonian possesses 
3 integrals of motion linear in the momenta (the three components of the angular 
momentum L) and a further 3 integrals of motion which are quadratic in the momenta 
(the components of the Laplace-Runge-Lenz vector A) : 

A = KxL + m^, (81) 

l x l 

Of course there are 2 independent relations among the seven constants of motion 
H, L, A, reducing to five functionally independent constants, which entitles the Kepler 
problem in S 3 to be maximally superintegrable. These relations are: 

2_„,2;„2 , fo™IT t2\t2 



A z = m A k z + [2mH - L z )L Z , A-L = (82) 

The algebra of Poisson brackets for the components of L and A has the following 
form: 

{Ai,Aj} = e ijk L k {2mH - 2L 2 ) , {U, A,-} = -e ijk A k . (83) 

Should these commutation relations be computed in a sphere with curvature k, 
the term — 2L 2 in the first Poisson bracket would appear as — 2kL 2 , making expressions 
dimensionally correct. This term displays clearly the effects due to the curvature of 
the configuration space, when it is compared (|83p with the commutation relations for 
the corresponding algebra in the flat case, which are 

{Ai,Aj} = eijkLk2mH, {Li, A,} = -e ijk A k . (84) 

From (f83|) we see that the components of L and A do not have the commutation 
relations of a Lie algebra, because of the presence of higher order terms in the Poisson 
bracket; the quadratic term L k H is already present in the flat case, but on the curved 
sphere S 3 cubic terms L k \? appear as well. In the flat Euclidean case, this algebra 
has been studied under the name of Higgs algebra [T2] . 

Then we consider this system as an example of the situation discussed in the 
previous section: the Kepler Hamiltonian admits two additional constant of motion 
which are also in involution: the components of the angular momentum anf of the 
Runge-Lenz vector along any fixed direction. If we take the z axis, then we get 
H, L%, A3 as three constants of motion in involution. 

Now we note that the Kepler problem has axial symmetry around any axis. In 
particular, the third constant A3 in the Kepler problem in Euclidean belong indeed 
to the family (J6j) with the values a = 0,(3 = 1,7 = for the parameters. While 
the coordinate system behind the family a ^ 0,(3 = 0,7 = were oblate spheroidal 
coordinates in E 3 and their analogous in S 3 , the family a ^ 0,f3 ^ 0,7 = turns out 
to be separable in parabolic coordinates in E 3 and in its analogous for S 3 . 

To end, we mention the possibility, which exists in S 3 , to rescale the Runge- 
Lenz vector by a factor which depends on the integrals of motion (so the property 
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of being a constant is not disturbed) in such a way that the Poisson brackets of new 
components, together with those of angular momentum closes a Lie algebra so(4). A 
rescaling ensuring this [35] is: 



R = A 



L 2 +m(\/H 2 + k 2 -H 



-1/2 



(85) 



(the rescaling factor is always positive, no matter neither the value nor the sign of H) 
and the components of L and R turn out to have the so(4) commutation relations 
with respect to the Poisson bracket (compare to ([55)) ): 

{Li, Lj} = —eijkLk , {Ri, Rj} = —tijkLk , {Li, Rj} — —CijkRk -(86) 

For these commutation relations (|5S|) the first Casimir reads R 2 + L 2 . Direct 
computations reveals that this Casimir is: 

R 2 + L 2 = m(H + ^H 2 + kA (87) 

Either by solving this equation for H or by a direct tedious but straightfoward 
computation, this allows us to write Hamiltonian (|80|) in terms of this Casimir: 



2mH = R 2 + L 



2 m 2 k 2 



R 2 



This discussion has been purely classical. We simply mention that this last 
property can be also discussed in the quantum case, where now the new global so(4) Lie 
algebra symmetry allows us to derive the energy spectrum first found by Schrodinger 
for the Kepler problem in S 3 [37J HH1 HS1 US] ■ Full details on this new derivation can 
be found in [25]; we must recall this approach is conceptually very different from the 
standard Pauli discussion for the Euclidean Kepler problem, where there is a different 
Lie algebra in each energy eigenspace, with its isomorphism class depending on the 
energy sign. 

5. Concluding remarks. 

We have discussed a special class of three dimensional completely integrable either in 
the flat Euclidean 3D space as well as in the standard sphere S 3 ; these systems are 
required to have a rotational symmetry axis, so they have axial symmetry. Further 
to the Hamiltonian H, there are two additional constants for such a system: L, the 
angular momentum with respect the symmetry axis, and Hi = l/(2rn)Pig % 2 ( x )Pj + 
>F(x). The general form for the tensor g^ix) which guarantees that the 'kinetic' parts 
of H and H2 commute among themselves and with L is given, and two special cases 
are mentioned; we discuss one of them, corresponding geometrically to separability 
in oblate spheroidal coordinates in full detail; the other corresponds to separation of 
variables in parabolic coordinates. Integrability determines the 'potentials' V(x), >F(x) 
to have a particular form in terms of undetermined functions of the coordinates, 
and the equations of motion are written as three equations each one involving one 
coordinate (separation of variables). In the quantum case, one of the three separated 
wave equations is trivial and each of the other two can be easily transformed into a 
spheroidal wave equation, for which the solutions have been studied. 

If no symmetry conditions are imposed a similar study can be performed. 
Separation of variables can be achieved in this case, as one could expect, in terms of the 
(general) ellipsoidal coordinates. In the quantum case, we obtain three similar wave 
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equations, one for each of the variables, each one in terms of a different function. These 
functions appear in the form of the potentials as a consequence of the integrability 
condition. 

Finally, we have carried the analysis of the systems with axial symmetry to the 
three dimensional sphere S 3 in the classical case. The results obtained are quite similar 
to the flat ones, in agreement to the idea that results in this area for the constant 
curvature spaces are essentially 'the same' as in the flat case. This is illustrated with 
the particular case of the three dimensional Kepler problem. Extension to higher 
dimensions should be also possible [5J [TU] 

Acknowledgements. 

Partial financial support is acknowledged to the Junta de Castilla y Leon Project 
VA013C05, the Ministry of Education and Science of Spain projects MTM2005-09183 
and FIS2005-03988 and Grant SAB2004-0169 and the Russian Science Foundation 
Grant 04-01-00352. G.P. was also supported by the programme ENTER-2004/04EP- 
48, E.U.-European Social Fund (75%) and Greek Ministry of development - GSRT 
(25%) and by RFFI grant 07-01-00234. 



References 

[l 

[2 



p: 

[6. 

[7. 

[»: 

[9 

[10 

[11 

[12; 

[13; 
[1< 

[is: 

[ie; 



A. Ballesteros, F.J. Herranz, M. Santander, T. Sanz-Gil, Maximal superintcgrability on 
TV-dimensional curved spaces, J. Phys. A 36, L93-99 (2003). 

S. Benenti, Intrinsic characterization of the variable separation in the Hamilton- Jacobi equation, 

J. Math. Phys. 38, 6578-6602 (1997). 
J.F. Carinena, M.F. Ranada, M. Santander, Central potentials on spaces of constant 

curvature: The Kepler problem on the two-dimensional sphere S 2 and the hyperbolic plane 

H 2 , J. Math. Phys. 46, 052702, 1-25 (2005). 
J.F. Carinena, M.F. Ranada, M. Santander, Superintcgrability on curved spaces, orbits and 

momentum hodographs: revisiting a classical result by Hamilton, J. Phys. A To appear (2007). 
M. Crampin, Conformal Killing tensors with vanishing torsion and the separation of variables 

in the Hamilton-Jacobi equation, Diff. Geom. Appl. 18, 87-102 (2003). 

B. Demircioglu, §. Kuru, M. Onder and A. Vercin, Two families of superintegrable and 
isospectral potentials in two dimensions, J. Math. Phys. 43, 2133 (2002). 

P. DOMBROWSKI. J. Zitterbarth, On the planetary motion in the 3-Dim standard spaces M% 
of constant curvature K, Demonstratio Mathematica 24, 375-458 (1991). 

N.W. Evans, Superintegrability in classical mechanics, Phys. Rev. A 41, 5666-76 (1990). 

M. Gadella, J. Negro and G.P. Pronko, J. Phys. A: Math. Theor., 40 (2007), 10791. 

M. Gadella, M. Ioffe, J. Negro, G.P. Pronko, Integrable systems in ellipsoidal coordinates, 
to appear. 

L. Garcia-Gutierrez and M. Santander, Levi-Civita regularization and geodesic flows for 
the 'curved' Kepler problem, arXlv: 0707. 3810 [math-ph] 

V.V. Gritsev and Yu.A. KurOCHKIN, The Higgs algebra and the Kepler problem in R3 J. 
Phys. A 33, 4073-4079, (2000). 

F.J. Herranz, A. Ballesteros, Superintegrability on Thrce-Dimcnsional Ricmannian 
and Relativistic Spaces of Constant Curvature, SIGMA (Symmetry, Integrability 
and Geometry: Methods and Applications 2, 010, 22p (2006). Available online at 
http : //www . ends . de/ j ournals/SIGMA/2006/ PaperOlO/ 

J.T. HORWOOD, R.G. McLenaghan, R.G. Smirnov, Invariant classification of orthogonally 
separable Hamiltonian systems in Euclidean space, Comm. Math. Phys. 259, (2005) 679- 
709. 

L. Infeld and A. Schild, A note on the Kepler problem in a space of constant negative 
curvature, Phys. Rev. 67, 121-122 (1945). 

E.G. Kalnins, G.C. Williams, W. Miller, G.S. Pogosyan, Superintegrability in three- 
dimensional Euclidean space, J. Math. Phys. 40, (1999) 708-725. 



Classical and quantum integrability in 3D systems. 15 

[17] E.G. Kalnins, J.M. Kress, W. Miller, G.S. Pogosyan, Complete sets of invariants for 

dynamical systems that admit a separation of variables, J. Math. Phys. 43, (2002) 708-725. 
[18] V.V. Kozlov AND A.O. Harin, Kepler's problem in constant curvature spaces, Celest. 

Mechanics 54, (1992) 393-399. 
[19] §. Kuru, A. Tegmen and A. Vercin, Intertwined isospectral potentials in an arbitrary 

dimension, J. Math. Phys. 42, (2001) 3344-3360. 
[20] A. A. Makarov, J. A. Smorodinsky, K. Valiev, and P. Winternitz, A systematic search for 

nonrelativistic systems with dynamical symmetries, Nuovo Cim. A 52, (1967) 1061-1084. 
[21] J. Meixner, F.W. Schafke, AND G. Wolf, Mathieu Functions and Spheroidal Functions and 

Their Mathematical Foundations, Springer- Verlag, 1980. 
[22] W. Miller, J. Patera and P. Winternitz, Subgroups of Lie groups and separatio of variables, 

J. Math. Phys. 22, (1981) 251-260. 
[23] L.M. NlETO, H.C. ROSU AND M. Santander, Hydrogen atom as an eigenvalue problem in 3D 

spaces of constant curvature and minimal length, Mod. Phys. Lett. A 14, 2463—2469, (1999). 
[24] M.N. Olevski, Triorthogonal systems in spaces of constant curvature in which the equation 

A3« + Xu = allows a complete separation of variables, Mat. Sb. 27 (69), 379—426 (1950) 

(In Russian). 
[25] G.P. Pronko, Kepler problem in the constant curvature space, arXiv:0705.3111 j2 [math-ph] 
[26] M. Santander and T. Sanz-Gil, From oscillator(s) and Kepler(s) potentials to general 

superintegrable systems in spaces of constant curvature, Rep. on Math. Phys., 55, pp. 371- 

383, (2005). 
[27] E. Schrodinger, A method of determining quantum mechanical eigenvalues and eigenfunctions, 

Proc. R.I.A. A 46, 9-16 (1940). 
[28] E. Schrodinger, Further studies on solving eigenvalue problems by factorization, Proc. Roy. 

Irish Acad. A 46, (1941) 183-206. 
[29] A.F. Stevenson, Note on the 'Kepler Problem' in a spherical space, and the factorization 

method of solving eigenvalue problems, Phys. Rev. 59, 842-843 (1941). 
[30] T.G. VOZMISCHEVA, Integrable problems of celestial mechanics in spaces of constant curvature, 

Astrophysics and Space Science Library, 295. (Kluwer Academic Pub., Dordrecht, 2003). 



